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CONVERGENCE AND OPTIMALITY OF HIGHER-ORDER 
ADAPTIVE FINITE ELEMENT METHODS FOR EIGENVALUE 

GLUSTERS 

ANDREA BONITO* AND ALAN DEMLOWt 

Abstract. Proofs of convergence of adaptive finite element methods for approximating eigen¬ 
values and eigenfunctions of linear elliptic problems have been given in a several recent papers. A 
key step in establishing such results for multiple and clustered eigenvalues was provided by Dai et. 
al. in [S], who proved convergence and optimality of AFEM for eigenvalues of multiplicity greater 
than one. There it was shown that a theoretical (non-computable) error estimator for which stan¬ 
dard convergence proofs apply is equivalent to a standard computable estimator on sufficiently fine 
grids. In Gallistl used a similar tool to prove that a standard adaptive FEM for controlling 

eigenvalue clusters for the Laplacian using continuous piecewise linear finite element spaces converges 
with optimal rate. When considering either higher-order finite element spaces or non-constant dif¬ 
fusion coefficients, however, the arguments of [8] and m do not yield equivalence of the practical 
and theoretical estimators for clustered eigenvalues. In this note we provide this missing key step, 
thus showing that standard adaptive FEM for clustered eigenvalues employing elements of arbitrary 
polynomial degree converge with optimal rate. We additionally establish that a key user-defined 
input parameter in the AFEM, the bulk marking parameter, may be chosen entirely independently 
of the properties of the target eigenvalue cluster. All of these results assume a fineness condition on 
the initial mesh in order to ensure that the nonlinearity is sufficiently resolved. 

Key words, eigenvalue problems, spectral computations, a posteriori error estimates, adaptiv¬ 
ity, optimality 

AMS subject classification. 65N12, 65N15, 65N25, 65N30 


1. Introduction. There has been high interest in recent years in the develop¬ 
ment and analysis of adaptive finite element methods (AFEM) for approximating 
eigenvalues and eigenfunctions of elliptic operators. We consider the following model 
eigenvalue problem: Find {uj,Xj) G x K. such that {uj,Uj) = 1 and 

a{uj,v) = Xj{uj,v), V £ iJQ(f2). (1.1) 

Here H C d = 2,3, is a polyhedral domain, a{u,v) := f^Vu ■ Vvdx and 
(u,v) := J^uvdx. There is then a sequence of eigenvalues 0 < Ai < A 2 < A 3 < .... 
and corresponding L 2 -orthonormal eigenfunctions iti,U 2 ,... satisfying (II. ip . Given a 
nested sequence of adaptively generated simplicial meshes {Ti}i>o with associated 
finite element spaces {Vf}i>o C the corresponding discrete eigenvalue 

problem is: Find Xij) £ x K such that (m£,j, Wj) = 1 and 

a{ue,j,v) = Xej{uej,v), v £ Yg. ( 1 . 2 ) 

We seek to approximate an eigenvalue cluster {Xj}j^j and associated invariant sub¬ 
space W := span{'Uj}jgj. Our index set J is given by J := {n-b 1, ...,n + N} for some 
n > 0, N > 1. The corresponding discrete sets are {Xij}j^j and Wg = spsLii{ugj}j^j. 
AFEM for eigenvalues are typically based on the standard loop 

solve —estimate -o mark —> refine. 
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To estimate the finite element error, AFEM employs local error indicators := 

SjeJ where r]g{T, X(j) is a standard residual error indicator for 

the residual —Auij — Let 0 < 0 < 1 be a given parameter, rj is used in mark 

to select a smallest set Ait cTe satisfying the Ddrfler (bulk) [10] criterion 

E ^ ^ E (1-3) 

TGMe TeTe 

Proofs of convergence and optimality of AFEM for approximating (IlB have been 
given in several papers. The first proof of optimality of AFEM for controlling simple 
eigenvalues and eigenfunctions was given in [9] . Other papers concerning convergence 
of AFEM for simple eigenvalues include [I11511I2]- The paper [5] contains a proof 
of optimal convergence of standard AFEM for an eigenvalue with multiplicity greater 
than one, while m proves a similar result for clustered eigenvalues. These papers 
mirror AFEM convergence theory for source problems (cf. |7]) in that they first prove 
that the AFEM contracts at each step. An optimal convergence rate dependent on 
membership of the eigenfunctions in suitable approximation classes is then obtained 
by standard methods. All require that the maximum mesh diameter in the initial 
mesh be sufficiently small to suitably resolve the nonlinearity of the problem. The 
behavior of AFEM for eigenvalues in the pre-asymptotic regime was studied in |13j . 
where the authors proved plain convergence results (with no rates) starting from 
any initial mesh. These results guarantee convergence of AFEM for general elliptic 
eigenproblems to some eigenpair, but not generically to the correct pair. 

The works of Dai et al. |8] and Gallistl m are most relevant to ours. In [Q the 
authors establish convergence of an AFEM for a multiple eigenvalue of a symmet¬ 
ric second-order linear elliptic operator for arbitrary-degree finite element spaces. A 
similar result is stated for eigenvalue clusters, but not all steps of the proof are pro¬ 
vided and the asymptotic nature of and constants in the results arising from the proof 
suggested in [8] depend on spectral resolution within the target cluster. Approxima¬ 
tion of eigenclusters of the Laplacian using piecewise linear elements is considered in 
m- The framework of m is cluster-robust, that is, all constants and the asymptotic 
nature of the estimates depending only on separation of the target cluster from the re¬ 
mainder of the spectrum. This leaves open the question of cluster-robust convergence 
results for AFEM for eigenvalue clusters using Lagrange spaces of arbitrary polyno¬ 
mial degree. We fill this gap by showing that standard AFEM for eigenvalue clusters 
using polynomials of arbitrary degree also converge optimally. While we consider only 
conforming simplicial meshes, we also provide a key step in extending such analysis to 
quadrilateral elements of any degree, meshes with hanging nodes, and discontinuous 
Galerkin methods; cf. |1] for analysis of the source problem. The analysis of m 
additionally does not immediately apply in the case of non-constant diffusion coef¬ 
ficients. In contrast, our results extend as in |8] to general symmetric second-order 
linear elliptic operators (see Remark [SG])- 

We briefly explain the difficulty which we resolve. A key step in standard AFEM 
convergence proofs establishing a certain continuity between error indicators on ad¬ 
jacent mesh levels. In the case of multiple or clustered eigenvalues, the ordering and 
alignment of the discrete eigenfunctions may change between mesh levels even on fine 
meshes. Thus {uij,Xtj) and (u£+ij, Xt+ij) may not approximate the same eigen¬ 
pair, making the comparison between r]t(T,utj, Xtj) and r/t^i{T,ui+ij, Xt+ij) used 
in standard AFEM convergence proofs irrelevant. A critical contribution was made 
in [5], where this problem was circumvented by first analyzing a theoretical (non- 
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computable) AFEM based on error indicators := f^j{T,Uj, Xj)^ aligned 

to the fixed continuous cluster, fxi may be viewed as an indicator for a “pseudo¬ 
residual” —AAiUj — XjP^Uj, where and Pt are projections onto Wf defined later. 
Given 9 G (0,1], the marking strategy consists of selecting the smallest set Ait C Pt 
of XI with 


> 9 Y, ^tt{T)\ (1.4) 

TGMe 

Optimality is guaranteed (see m) provided the initial mesh is fine enough and 9 < 9* 
for some 9* < 1 only depending on Xl and the initial subdivision To ■ To deduce optimal 
properties of the practical algorithm based on r]t(T) instead of ^t{T), an elementwise 
equivalence of rjt and ^t is used to show that the corresponding AFEM are essentially 
equivalent. The extension of the proofs given for this final step in Hill] to the approx¬ 
imation of clustered eigenvalues using higher-degree finite element spaces is however 
unclear. We provide (Lemma 13.111 a more direct proof of the equivalence between 
pit and rjt which is valid for finite element spaces of arbitrary degree. Extension of 
the AFEM convergence results for eigenvalue clusters of m to higher-degree finite 
element spaces is then immediate (see Theorem l3.5l and Corollarv l3.9l) . Our estimator 
equivalence argument could also be used with minimal modification to complete the 
proof of AEEM optimality for eigenvalue clusters given in [8] Theorem 7.1]. It is 
however not immediately clear how to modify certain other estimates in [8] in order 
to obtain cluster-robust results. 

Our approach also yields equivalence constants between jit and rjt which are 
independent of essential quantities on sufficiently refined meshes, leading to a final 
result which is more robust with respect to the properties of the target cluster. In view 
of (|1.3I1 and (|1.4I1 . optimal AEEM require 9 <9* := C9*, where C < I is independent 
of N, A„+i, and Xn+N (see Remark [3.81) . In contrast, the analysis of [TT] implies 
that a stricter (smaller) choice of 9 may be necessary as N and the ratio A„+Ar/A„+i 
increase. As numerical experiments confirm (]J4]), this improvement is important from 
a practical standpoint as it establishes that no knowledge of cluster properties is 
needed in order to choose 9 correctly. It also potentially increases computational 
efficiency by confirming that 9 may be chosen reasonably large even when computing 
large clusters. On the other hand, m implies that 9 should be taken to be small as N 
increases, thus yielding a large number of adaptive loops and decreasing computational 
efficiency. 

The remainder of the note is structured as follows. In 1}2] we give further as¬ 
sumptions and preliminaries. In ]J3] we prove our main result and briefly sketch its 
application in the proof of AEEM convergence for eigenvalue clusters. A brief set of 
numerical experiments in ]2] illustrates our theoretical results. Finally note that we 
do not give a full proof of eigenvalue convergence, as most of the analysis of m holds 
verbatim for higher-order elements given the equivalence of the theoretical and prac¬ 
tical error estimators. We also largely employ the notation of that work. Familiarity 
with m is thus essential to the reader. 

2. Preliminaries. 

2.1. Finite element meshes and spaces. Let 17 be a polyhedral domain and 
denote by T the set of all conforming refinements of To obtainable by the newest- 
vertex bisection algorithm typically used in AFEM or its generalization to higher 
dimensions; cf. [1II711I11IS]. By construction all T £ T are uniformly shape regular. 
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Given T € T, we denote by V(T) C Hq{Q) the space of continuous piecewise 
polynomials of arbitrary but fixed degree r subordinate to T. For each £ > 0, we 
abbreviate Vf := Y{Ti). The AFEM produce sequences of subdivisions {Te}e>o C T 
such that Ti+i is a refinement of Ti- Thus Ye C V^+i. 

We also define several projection operators. Let Pe : L 2 {^) —>■ be the L 2 

projection onto W^, and let Gi : Hq{Q) Ye be the Ritz projection. That is, 

{Peu,v) = (u,v), v€We, a(Geu,v) = a(u,v), v €Ye- ( 2 . 1 ) 

In addition, we define Ae := PeoGe- Finally, P : ^ 2 ( 1 !) —>■ IF is the L 2 projection onto 
the continuous invariant subspace W corresponding to the target eigenvalue cluster. 

2.2. Error estimators and AFEM. Given a discrete eigenpair {uej, Xej) with 
j € J and T G Te-, the computable local error indicator r]e(T) is given by 

miT)^ ■= (/lyll+ Afj-Mfjj|^ +j-Jllgji) , TGTe- (2-2) 

j&J 

Here !•] is the jump across the element boundary, || • ||£) is the L 2 norm over D 
and hx ■= diam(T). Let ||| - ||| = y/a{-, •) be the energy norm. The computable 
error estimator := reliably controls ~ up to 

higher-order terms which can be proved to be negligible under the assumption that 
Hq := maxj’gT^ hx is sufficiently small. 

The theoretical local indicator fie{T) is given by 


^leiTf := ^ {h^X.PeUj + AA^n,||2 + hxUYAeUjjWlx) ■ (2.3) 

j&J 

A slight modification of Proposition 4.1 of m shows that J^xeTe /^(^)^ bounds 
|||it — A^tt||p up to higher-order terms which are negligible if Hq is sufficiently small. 
In m the volume term is simplified to h^\\XjPeUj\\^ due to the use of piecewise 
linear finite element spaces, and the definition in the higher-order case is left open 
with our choice given as one possibility (cf. Remark 9 (f) of [H]). Our definition of ne 
also differs from that given in [8]. The difference does not appear to be essential for 
purposes of proving estimator equivalence. However, certain other estimates needed 
to obtain a cluster-robust final result have thus far only been obtained using the choice 
of fie suggested in m and used here. 

2.3. Assumptions on the initial mesh. Proofs of optimal convergence rates 
for adaptive FEM require assumptions on the eigenspaces and their resolution by the 
initial mesh. First, the continuous eigencluster must be separated from the remainder 
of the spectrum, i.e., A„ < A„+i < ... < A„+ 7 v < A„+Ar-i-i (here we take Aq = 0). 
In addition, the discrete cluster respects this separation in that Xe^n < A„+i and 
Xe,n+N < A„+JV+ 1 . Let then AIj := sup^-^gT maXjg{i....,dimV 4 \J maxfcgj 
The implied requirement on the initial mesh may be more or less strict depending on 
the separation of the cluster from the remainder of the spectrum. 

Reliability of our a posteriori estimators is only guaranteed for sufficiently fine 
initial mesh Tq. We have for example from Proposition 4.1 of m that 

fieiTf + A|(l + MjfH^Q^lWu, - Aeu,\\\^ 

-TeTf 


Uj — AfUjlll^ < C 


(2.4) 
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Thus the estimator is reliable when Hq is small enough to guarantee that CXj{\ + 
Mj)‘^Hq‘’ < This assumption is a priori and cannot be rigorously verihed com¬ 
putationally. Any assumption on the resolution of eigenvalues within the cluster is 
however avoided. Cf. [niiiiiis] for more discussion. 

3. Main result. 

3.1. Equivalence of theoretical and practical estimators. Below is our 
main result. 

Lemma 3.1 (Estimator Equivalence). Assume that Hq = maxTeTj hx is small 
enough so that 

max \\uj — KiUjWo, < \/l + {2N)-^ — 1, j £ J. (3.1) 

jeJ 

Then for T 

< ^77,(r)2 < 3A*,(r)2. (3.2) 


Proof Using (11.11) . (11.21) . and (EH), we compute 

m^J m^J 

mGJ m^J 

— ^ ^ ^£,rn{^£Uj ■) U£,rn)ui^rn’ 
mGJ 

Also, 

m£j 

Combining the above two equations yields 

m£j 

[ va , m,i = ^ {Aiuj,ui 

,m )IVm^ ,m| • 

m^J 


(3.3) 


(3.4) 


(3.5) 


We now define the vectors V := [XjPiUj + AA^Uj], Vi := [Xe^mUt^m + Aui^m], 
W := [iVAfMj]], and Wi := [[[Vm^j]]. Here n -b 1 < j, m < n -|- iV. These quantities 
vary with x £ Tl (for V and Vi) or x lying on the mesh skeleton (for W and Wi), 
but this fact has little significance for the time being. We also define the matrix 
M := [{AiUj,Ui^rn)]i with j the row and m the column index. The relationships (13.5|) 
may now be written 


V = MVi, W = MWi. (3.6) 

Let II • II 2 denote the operator £ 2 (R^)-norm and | • | the Euclidean length. The proof 
of Lemma 13.11 is thus reduced to showing that 

||M||2 <vW2 , \\M-^\\2<V2 


(3.7) 
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under the stated hypotheses. 

We now employ several basic facts from linear algebra. Letting v G , we have 
that IlMulll = v^M*Mv < ||M*M|| 2 |up and thus ||M ||2 < ||M*M|| 2 . Because M*M 
is symmetric and positive, ||M*M|l 2 is equal to the maximum eigenvalue of M*M. 

A short computation shows that B := MM* = A^rim)]. A standard 

result from linear algebra is that AC and CA are isospectral for square A, C, so B is 
isospectral with M*M. Therefore we have that ||M*M ||2 = ||B|| 2 , and both quantities 
are equal to the maximum eigenvalue of B. The matrix B was analyzed extensively 
in the proof of Lemma 5.1 of m- In particular, B is nonsingular under the condition 
(EU, and by (5.2) and following of [11], we have 


2iV - 1 2iV + 1 

< B, ,, < 


2N 


2N 


and 




N-1 
2N ■ 


Gershgorin’s theorem thus gives that the eigenvalues {ai} oi B satisfy 

1 2A^- 1 

2 


2N 


N-1 2N+1 

< CFi < 


2N 


2N 


N-1 3 , , 

-= -, 1 <i< N. 

2N 2’ - - 


(3.8) 


(3.9) 


Thus |1M||| < ||M*M ||2 = ||MAf *||2 < |, which is the first inequality in (13.71) . 

The invertibility of B guarantees the invertibility of M. Computing as above 
yields ||M “^||2 < ||(M“^)*M “^||2 = ||i?“^|l 2 . Because B is positive and diagonaliz- 
able, we have from (13.91) that ||i?“^|| = mini<^<ff o- — completing the proof of 

the second inequality in (Irfl) . □ 

Remark 3.2 (Mesh Fineness Assumption). In it is shown for piecewise 
linear elements that under the assumption m, 


ti^T? < N{K+M/Xn + < (A„+m/A„+i)"(21V2 ^ 4 ^r 3 )^^(T) 2 . (3.10) 


The first bound in (13.101) holds with no restriction on Hq . Our proof thus yields an 
improved cluster-independent bound under the fineness assumption (lO) . (lO) also 
gives an improved bound for r]i in terms of . Significantly loosening the restriction 
(ED) appears to be substantially more difficult. M is invertible for e < Vl N ^ — 1, 
but the obtained upper bound for ||M “^||2 degenerates as e f v*! + N~^ — 1. The 
condition e < Vl + N~^ — 1 is not substantially weaker than EH), and control over 
the equivalency constants is lost as e f Vl + N~^ — 1, so we retain EH from m- 
Remark 3.3 (Alternate proof of estimator equivalence). In Lemma 3.3 of it 
is shown that if An+i = Xn+N and Yi is of arbitrary polynomial degree, 


fi{Tf < Npi{Tf < N^C p{Tf + EE \Xij — A^^fcP^TlI^^.fellT 1 ) (3-11) 

\ jejkej j 

where C —>■ 1 as Hq —>■ 0 and fii is similar to In our analysis the obtained bounds 
are independent of N and the double sum term on the right is completely absent. The 
absence of the double sum term on the right is important, as it is of higher order for 
multiple but not for clustered eigenvalues; cf. Remark 9 (f) of Ulf . Usable results 
for clusters may be obtained by summing (13.111) over each multiple eigenvalue in a 
cluster, but the constants in and asymptotic nature of the resulting estimate then are 
not cluster-robust. 

Remark 3.4 (General symmetric operators). In AFEM optimality is ana¬ 
lyzed for eigenvalues of arbitrary second-order symmetric operators —V • (AVu) + cu. 
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Letting —+ our equiv¬ 

alence analysis extends immediately to the general ease. We consider only A = I, 
c = 0 because analyzing optimality in the general case involves consideration of data 
oscillation, which would make our presentation much more involved; cf 

3.2. AFEM optimality. Let 


■= {v = (vi,...,vj) e HKVlY : |vUj < oo} , 
where Q < a < r/d and 


IvUj := sup m’^ inf inf 

mSN T£T,#T—#7o<m 




1/2 


Given Lemma |3.11 the analysis of m mostly holds verbatim for higher-order finite 
element spaces. In particular, the following counterpart to Theorem 3.1 of m holds. 

Theorem 3.5 (Quasi-Optimal Approximation of Eigenspaces). Assume that the 
bulk parameter 6 < 1 and initial mesh size Hq are sufficiently small. Assume that for 
some 0 < cr < r/d, u := {un+i, ...,Un+N) S Af,. Then there exists C > 0 depending 
possibly on To but independent of other essential quantities such that 

fg|||n,-A,n,|||2j < C'(#T£ - #ro)-nuU- (3.12) 


Several remarks are in order. 

Remark 3.6 (Approximation classes for clusters). We consider approximation 
classes Af, for the entire cluster of eigenfunctions, thereby imposing that all the eigen¬ 
functions in the cluster can be approximated simultaneously (on the same subdivision) 
with rate a. It is equivalent to requiring that each eigenfunction belongs to 

Act :=-1G iLn(n) : RU., := sup inf inf III?; — ())||| < oo 

i uv ; I I reT,#r-#ro<m,/ev(r) 

but the norm equivalence constant depends on J. Our choice of approximation class 
guarantees that the constant in estimate (13.121) is independent of J. 

Remark 3.7 (Equivalent approximation classes). In \llf optimality is expressed 
with respect to the class 

Ba:= IvG : \v\b„ ■= sup inf || V?; - PD(r) V?;|| < oo 

( mGN TeT,#T-#To<m 

where Pd(T) the L 2 -projection onto D(T), the space of piecewise polynomials of 
degree r — 1 subordinate to T. Employing this class requires proving equivalence be¬ 
tween Act (see E.emark \S.6\} and Ba, that is, approximation by functions in ©(T) is 
equivalent to approximation by gradients of functions in V(T). In \lTf . the needed 
result is obtained by citing in which the equivalence is shown up to data oscil¬ 

lation terms. However, the necessary equivalence has recently been shown to hold on 
arbitrary meshes without data oscillation terms in m, which could simplify the proof 
of Proposition 3.1 of m- 

Remark 3.8 (Robustness with respect to the cluster size). The dependence 
of Hq and 0 on various quantities is given in detail in m- The threshold condi¬ 
tion for 9 depends on the second equivalence constant in (|3.2I) . Thus m requires 
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0 <0 < [C'2(A„+Ar/A„+i)^(2A^^+ 4A^^)] ^ (cf. Lemma 7.3), where C 2 is independent 
of essential quantities (especially cluster properties). We only require 0 < 6 < 
thus establishing that rate-optimal adaptive convergence holds with the marking pa¬ 
rameter chosen independent of the cluster properties. The analysis of m however 
still indicates that the asymptotic nature of the optimality result may be more or less 
pronounced depending on cluster properties (cf. (12.4p ). 

Finally, optimal approximation of eigenvalues follows as in m Corollary 3.1]. 

Corollary 3.9 (Quasi-Optimal Approximation of Eigenvalues clusters). Under 
the conditions of Theorem \3.5\ there exists Cn,N > 0 such that for n-\-l<j<n-\-N, 

\Xe,j -Xj\< Cn,N{#T, - #ro)-^" ^ |u,|l. 

jeJ 


Note that the constant Cn,N above depends on properties of the eigenvalue cluster; 
cf. m Corollary 3.1] for a precise expression. 

4. Numerical experiments. To illustrate our result, consider the slit domain 


O :=(-!, 1)2 \ 


f conv{(0.5,0), (1,0)} U conv{(0,0.5), (0,1)}U 
yconv{(—0.5, 0), (—1,0)} U conv{(0, —0.5), (0, —1)} 


subdivided into 16 squares, each of diameter 0.5. The finite element spaces con¬ 
sist of continuous piecewise polynomials of degree r = 1, 2, or 3 in each coordi¬ 
nate direction. The adaptive algorithm is implemented within the deal.ii [2] library. 
We take n = 0 and consider the cases N = A, 12. Our highest-precision compu¬ 
tations indicate that the first twelve eigenvalues (truncated to 6 decimal places) 
are 10.147392, 17.662596, 17.662596, 19.739208, 26.101811, 38.349159, 38.349159, 

46.553966, 49.149607, 49.149607, 49.348022, and 49.348022. In FigureliH we plot the 

1 /2 

error estimators pi := against the total number of degrees of freedom 

for different values of marking parameters 0 and cluster sizes N. 

We do not know the exact eigenfunctions and so cannot plot actual errors. How¬ 
ever, in the case of the Laplacian m guarantees plain convergence (without rates) 
of AFEM for simple and multiple eigenvalues to the corresponding continuous eigen- 
pairs. These results rely on completely differently proof techniques than do ours, and 
adaptation to the case of clustered eigenvalues is straightforward. The analysis of m 
also guarantees that /it —t 0 as £ —>■ 00 , which yields reliability of pi for i 

sufficiently large (cf. (12.41) 1. It is thus meaningful to track p£ instead of actual errors. 
In addition, employing a cluster-independent marking parameter as suggested by our 
theory is reasonable also in the pre-asymptotic range as the plain convergence analysis 
of |13j guarantees that the algorithm will eventually reach the asymptotic range. 

Optimal decay rates of —r/2 are observed provided 0 < .95 when r = 2,3 and 
0 < 0.9 when r = 1. According to Corollary 13.91 (cr = 1/2), this indicates an optimal 
decay rate of — r for the error in approximating each eigenvalue in the cluster. As 
guaranteed by our analysis, the range of 0 for which optimality is recovered is not 
affected by the cluster properties. In contrast, following m as explained in Remark 
I3.8l leads to the restriction 0 < [C'2(A„+Ar/A„-|-i)^(2N2-|-4N^)]“^ with C 2 independent 
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Fig. 4.1. Values of rji versus dim(Yi) during the adaptive process for marking parameters 
0 = 0.8, 0.85, 0.9, 0.95, 1, with n = 0, N = A (top row) and n = 0, N = 12 (bottom row) in each cases 
using continuous piecewise polynomials of degree r (column r). The optimal decay rate of —r/2 is 
observed as soon as 0 < 1 when r = 2, 3 and 9 < 0.95 when r = 1. 


of essential quantities. For our particular computations, this yields: 


e < cV 
e < 


19.739208\ 


-4 


10.1473927 2-42+ 4 -43 


2.42 X 10“^a 


— 4/^—1 


= 0, iV = 4, 


-4 


_ ,.2.48X10-C-, n = 0,iV = 12. 

10.1473927 2-122 + 4-123 2 > 


Thus following precisely the theory of m would lead to a thousand-fold reduction 
in 6 when moving from our first to our second computational example. This would 
potentially require a massive increase in the number of AFEM iterations required 
in order to achieve a given error reduction. We have demonstrated theoretically and 
conhrmed computationally that this increase in computational expense is unnecessary. 
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